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We discuss the asymptotic dynamical evolution of spatially homogeneous brane-world cosmolog- 
ical models close to the initial singularity. We find that generically the cosmological singularity 
is isotropic in Bianchi type IX brane-world models and consequently these models do not exhibit 
Mixmaster or chaotic-like behaviour close to the initial singularity. We argue that this is typical 
of more general cosmological models in the brane-world scenario. In particular, we show that an 
isotropic singularity is a past-attractor in all orthogonal Bianchi models and is a local past-attractor 
in a class of inhomogeneous brane-world models. 
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^ . I. INTRODUCTION 

D 

' Higher-dimensional gravity theories inspired by string theory, in which the matter fields are confined to a 3- 
dimensional 'brane-world' embedded in 1 -I- 3 -I- d dimensions while the gravitational field can also propagate in the d 
extra dimensions (i.e., in the 'bulk') Q, are currently of great interest. The d extra dimensions need not be small or 
even compact in these theories, which effectively become 4-dimensional at lower energies. In recent work Randall and 
Sundrum Q have shown that for d — 1, gravity can be localized on a single 3-brane even when the fifth dimension is 
^ infinite. In this paradigm, Einstein's theory of general relativity (GR) must be modified at high energies (i.e., at earl 
times). An elegant geometric formulation and generalization of the Randall-Sundrum scenario has been given 
Much effort is has already been devoted to understand the cosmology of such a Randall-Sundrum-type brane-world 
scenario l^],^. 

The dynamical equations on the 3-brane differ from the GR equations by terms that carry the effects of imbed- 
ding and of the free gravitational field in the five-dimensional bulk The local (quadratic) energy- momentum 
. corrections are significant only at very high energies. In addition to the matter fields corrections, there are nonlocal 
' effects from the free gravitational field in the bulk, transmitted via the projection of the bulk Weyl tensor, that 
, contribute further corrections to the Einstein equations (e.g., the Friedmann equation) on the brane. Due to its 
O I symmetry properties, £fj_i, can be irreducibly decomposed (with respect to a timelike congruence u'') into 
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£t.u = -[-) [U {uf^u^ + \h^^) + Vf,^ + 2Q(^u^)] , (1) 



in terms of an effective nonlocal energy density on the brane, U, arising from the free gravitational field in the bulk, 
an effective nonlocal anisotropic stress on the brane, T^^jy, and an effective nonlocal energy flux on the brane, 

In general, the conservation equations do not determine all of the independent components of on the brane. 
In particular, there is no evolution equation for V^i, and hence, in general, the projection of the 5-dimensional field 
equations onto the brane does not lead to a closed system. However, in the cosmological context, in which the 
background metric is spatially homogeneous and isotropic, we have that 

- Q,, - P^. - , (2) 

where is the totally projected part of the brane covariant derivative. Since V^u = 0, in this case the evolution 
of £fj_i, is fully determined |2j . In general U — U{t) ^0 (and can be negative) in the Friedmann background . 
For a spatially homogeneous and isotropic model on the brane equation (|^) follows, and similar conditions apply 
self-consistently in other Bianchi models The Friedmann brane-world models have been extensively investigated 
il- 

There are many reasons to consider the classical dynamical evolution in more general spatially homogeneous Randall- 
Sundrum-type brane-world cosmological models, particularly in the early Universe in which the dynamical behaviour 
can be completely different to that of GR. In this article we shall consider the local dynamical behaviour of the general 
Bianchi type IX models close to the singularity. 

Due to the existence of monotone functions, it is known that there are no periodic or recurrent orbits in orthogonal 
spatially homogeneous Bianchi type IX models in GR. In particular, there are no sources or sinks and generically 
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Bianchi type IX models have an oscillatory behaviour with chaotic-like characteristics, with the matter density becom- 
ing dynamically negligible as one follows the evolution into the past towards the initial singularity. Using qualitative 
techniques, Ma and Wainwright (see 0) have shown that the orbits of the associated cosmological dynamical system 
are negatively asymptotic to a lower two-dimensional attractor. This is the union of three ellipsoids in consisting 
of the Kasner ring joined by Taub separatrices; the orbits spend most of the time near the self-similar Kasner vacuum 
equilibrium points. More rigorous global results are possible. Ringstrom has proven that a curvature invariant is 
unbounded in the incomplete directions of inextendible null geodesies for generic vacuum Bianchi models, and has 
rigorously shown that the Mixmaster attractor is the past attractor of Bianchi type IX models with an orthogonal 
perfect fluid 

All spatially homogeneous models in GR expand indefinitely except for the Bianchi type IX models. Bianchi type 
IX models obey the "closed universe recollapse" conjecture whereby initially expanding models enter a contracting 
phase and recollapse to a future "Big Crunch" . All orbits in the Bianchi IX invariant sets are positively departing; in 
order to analyse the future asymptotic states of such models it is necessary to compactify phase-space. The description 
of these models in terms of conventional Hubble- or expansion-normalized variables is only valid up to the point of 
maximum expansion (where H = 0). An appropriate set of alternative normalised variables, which leads to the 
compactification of Bianchi IX state space, were suggested in Q (see section 8.5.2) and have been utilized in the 
qualitative study of locally rotationally symmetric perfect fluid models [|l^ , of scalar field models with an exponential 
potential and scalar field models with a perfect fluid ||l2[. In particular, spatially homogeneous Bianchi models 
containing a scalar field (p with an exponential potential of the form V{(j)) = Aexpfc^ were studied in fl^ . For Bianchi 
type IX models, if the parameter k < \f2 then there exists a set of ever-expanding anisotropic Bianchi type IX models 
that isotropize and inflate towards an expanding power-law inflationary solution, except for the subset that recollapse 
I p^ , and if fc > \f2 then initially expanding Bianchi type IX models do not isotropize towards an ever-expanding 
isotropic model |]lT| . 



II. GOVERNING EQUATIONS 



The field equations induced on the brane are derived via an elegant geometric approach by Shiromizu ct al. 
using the Gauss-Codazzi equations, matching conditions and Z2 symmetry. The result is a modification of the standard 
Einstein equations, with the new terms carrying bulk effects onto the brane: 
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The brane energy-momentum tensor for a perfect fiuid or a minimally-coupled scalar field is given by 



(3) 



(4) 



(5) 



where u'^ is the 4-velocity, p and p are the energy density and isotropic pressure, and the projection tensor h^^^ = 
g^L' + u^u^ projects orthogonal to u^. We shall assume that the matter content is equivalent to that of a non-tilting 
perfect fiuid with a linear barotropic equation of state for the fluid, i.e., p = (7 — l)p, where the energy conditions 
impose the restriction p > 0, and the constant 7 satisfies 7 G [0, 2] from causality requirements. For a minimally 
coupled scalar field the energy density and pressure are, respectively. 



(6) 



and the conservation law is equivalent to the Klein-Gordon equation. A dynamical analysis of scalar field models 
indicates that at early times the scalar field is effectively massless 12 1. A massless scalar field is equivalent to a perfect 
fluid with a stiff equation of state parameter 7 = 2. In the early Universe, and close to the singularity, we expect that 
7 > 4/3. 

The bulk corrections to the Einstein equations on the brane are of two forms: flrstly, the matter flelds contribute 
local quadratic energy-momentum corrections via the tensor 5*^1/, and secondly, there are nonlocal effects from the 
free gravitational fleld in the bulk, transmitted via the projection S^^, of the bulk Weyl tensor. The matter corrections 
are given by 
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which is equivalent to 



Si_iu = YjP^UpU^ + -^p{p + 2p) hf_,^ 



(7) 



for a perfect fluid or minimally-coupled scalar field. The quadratic energy-momentum corrections to standard GR 
will be significant for K^p^ K^p in the high-energy regime. 

All of the bulk corrections may be consolidated into an effective total energy density, pressure, anisotropic stress 
and energy flux, as follows. The modified Einstein equations take the standard Einstein form with a redefined 
energy-momentum tensor: 
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where 



^2 ^A""^ ^2 



is the redefined perfect fluid (or minimally coupled scalar field) energy-momentum tensor with 
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where we recall that D^U = = V^i, = in the cosmological case of interest here. 

As a consequence of the form of the bulk energy-momentum tensor and of Z2 symmetry, it follows 
brane energy-momentum tensor separately satisfies the conservation equations, i.e., 

V^T^, = , 

whence the Bianchi identities on the brane imply that the projected Weyl tensor obeys the constraint 



(8) 
(9) 

(10) 
(11) 

that the 
(12) 
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In the spatially homogeneous case, the conservation equations for a non-tilting perfect fluid with a linear barotropic 
equation of state (or massless scalar field) reduce to 



p + 3H{p + p)^Q, 



(14) 



where the Hubble parameter H = d/a gives the volume expansion rate, and the nonlocal conservation equations 
reduce to an evolution equation for U : 



U + AHU = . 



(15) 



In the Friedmann case Eq. ( p^ yields the 'dark radiation' solution U =Uo {no/a) ■ 

The field equations are given by Eqns (|^) and (||). The generalized Friedmann equation, which determines the 
expansion of the universe or the Hubble function, in the case of spatially homogeneous cosmological models is 
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= -n^p ( 1 + -£-)- --^i? + + Ia + ^ 



2U 



(16) 



where is the scalar curvature of the hypersurfaces orthogonal to the fluid flow, which we associate with the 
cosmological fluid, and 2a^ = cr^^dab is the shear scalar in terms of the shear a^^. 

For a perfect fluid or minimally-coupled scalar fleld, the generalized Raychaudhuri equation on the brane becomes [El: 



H + H^ + lcj^ + ln\p + ip) [n^p{2p + 3p) + ITU] 



(17) 



where we have assumed (H), and in particular V^j,i, — 0, self-consistcntly on the brane 
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III. BIANCHI IX MODELS 



Following 0,|T^ we define Hubble- normalized (and dimensionless) shear variables S]+,I]_, curvature variables 
Ni, N2, N3 and matter variables ^i,Pi (relative to a group- invariant orthonormal frame), and a logarithmic (di- 
mensionless) time variable, r, defined by dr = Hdt. These variables do not lead to a global compact phase space, but 
they are bounded close to the singularity . The governing evolution equations for these quantities are then 

- (-7 - 2)E+ - 5+ (18) 
E'_ = (q - 2)S_ - (19) 
N[ = {q- 4I]+)iVi (20) 

Af^ = ((7 + 2E+ + 2V3S_)iV2 (21) 
iV^ = (g + 2S+-2\/3I]_)7V3, (22) 

where a prime denotes differentiation with respect to r, and K, S+ and S- are curvature terms that are defined as 
follows: 

K=^ {{Nl + Nl + N^) - 2{mN2 + N2N3 + N1N3)) (23a) 

S+ = l {N2 - N^f _ 1 TVi (2 TVi - 7V2 - N3) (23b) 
6 6 

S- = lV3{N2- N3) i-m +N2+ N3) (23c) 
6 

The quantity q is the deceleration parameter given by 

(Z EE 21]^ -I- 2S2_ + 1^0, + ^^ P,. (24) 
The decoupled Raychaudhuri equation becomes 

H' = -{l + q)H. (25) 
In addition, the generalized Friedmann equation reduces to the constraint 

+ K. (26) 

Due to the symmetries in the dynamical system, we can restrict ourselves to the set A^i > 0, -/V2 > 0, and N3 > 0, 
without loss of generality. (Note, for Bianchi IX models all of the iV^'s must be of the same sign.) We again note that 
these new normalized variables, for which the evolution equation for H has decoupled from the remaining equations, 
are bounded close to the initial singularity. 

In the above 

^2 ^tot ,^2 ptot 

where the fti (2 = 1 — 4) are given by 

2 

^1^^ = ^; P = h-m (28a) 



3iJ2' 



2 2 

n2 = nb = Pb = (27 - l)nb (28b) 



2U 1 

^3 = - ^^2^; Pu - flu (28c) 
A 



which satisfy the equations 



= = — ; Pa = -r!A, (28d) 



Vl[^[2{l + q)-iV,]U,, (29) 
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where 

Ti = 7,r2 = Tb = 27,r3 ^Tu = ^,Ti = TA = 0. (30) 

A minimally coupled homogeneous scalar field = with an exponential potential can also be included, where 
the energy-momentum tensor is given by 



TJ^ = cl^-.aCl^-.b - gab ( ^-.c^'' + Viq^) ] . (31) 
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The usual (linear) terms can be included via Q 

r!^ = ^2^$,P^ = ^2_^^ (32) 

where 

The evolution of these variables is then given by the conservation or Klein-Gordon equation, which can be written as 

= (g-2)^'- ^fc$2 (34) 

= $ ^(l + g) + ^fc*^ (35) 

The local quadratic corrections flb^ and Pb^ can be included via <i>b oc and ^Pf, cx and a generalized (quadratic) 
Klein-Gordon equation. Since an extensive analysis of scalar field models has shown that close to the initial singularity 
the scalar field must be massless ||l^ , |l3|] (see also |l5| , |l6| ), it is plausible that scalar field models can be approximated 
by a stiff perfect fluid close to the initial singularity (particularly regarding questions of stability) . 



A. Initial Singularity 

If > and -I- 3p*°* > (which follows from the various energy conditions) for all i < t^, then from the 
generalized Raychaudhuri equation (|l^) (or (|2^)) and using the generalized Friedmann equation (|l^) (or (^6|)) and 
the conservation equations, it follows that for a'o > (where ag = a(io)) there exists a time tf, with tb < to such that 
aitb) = 0, and there exists a singularity at tb, where we can rescale time so that tb ~ and the singularity occurs at 
the origin. 



From Eqn (17) we can find the precise constraints on A and Ua in terms of oq in order for these conditions to be 
satisfied at t — to- It then follows from the generalized Raychaudhuri equation, the generalized Friedmann equation 
and the conservation equations that if these conditions are satisfied at t — Iq they are satisfied for all < t < io: ^^id 
a singularity necessarily results. These conditions are indeed satisfied for regular matter undergoing thermal collapse 
in which the the local energy density and pressure satisfy p{2p + 3p) > (and is certainly satisfied for perfect fluid 
matter satisfying the weak energy condition p > and a linear barotropic equation of state with 7 > 1). On the other 
hand, it is known that a large positive cosmological constant A or a significant negative nonlocal term U counteracts 
gravitational collapse and can lead to the singularity being avoided in exceptional circumstances. 

Since the variables are bounded close to the singularity, it follows from Eqn (|2j) that < q < 2. Hence, from Eqn 
(^), H diverges as the initial singularity is approached. At an equilibrium point q~q*, where q* is a constant with 
< g* < 2, so that from Eqn (|2|) we have that H ^ {I + q*)~^t-^ as t ^ (r ^ -00) 



From Eqns (pj), \QM and the conservation laws it then follows that p ^ cx) as t — > 0+. It then follows directly 



from the conservation laws (29) that fib dominates as i ^ 0+ and that all of the other Vli are negligable dynamically 
as the singularity is approached. The fact that the effective equation of state at high densities become ultra stiff, so 
that the matter can dominate the shear dynamically, is a unique feature of brane cosmology. 
Note that the function defined in Bianchi IX models by 

Z={N^N2N^)\ (36) 



5 



satisfies the evolution equation 

Z' = 6gZ, (37) 

and is consequently monotone close to the singularity 0. 

In summary, generically the Bianchi type IX models have a cosmological initial singularity in which p — s- oo, and 
consequently fi;, dominates, as t ^ 0+. This can be proven by more rigorous methods jl^,^. It can also be shown 
by qualitative methods that the spatial 3-curvature is negligable at the initial singularity and from a comprehensive 
analysis that at later times p decreases and the low density regime in which GR is valid ensues. 



B. The Isotropic Equilibrium Point 

Hence we have shown that close to the singularity = p^, so that X) — ^f> ^^"^ ~ where 



(27 — Consequently, Eqn (26) can be written as 

VLy, = \- Y.\-Y?_ - (38) 

which can be used to eliminate fib from the governing equations. In particular, Eqn ( |2^ ) becomes 

g = 3(1 - 7)(S^ + S'_) + (37 - 1)(1 - K), (39) 

and the governing equations are given by the dynamical system (18-22), where q now given by the above expression. 

There is an equilibrium point of the dynamical system, denoted by given by E4. ~ S]_ = 0, and A^i = iV2 = 
A'3 = 0, which corresponds to spatially homogeneous and isotropic non-general- relativistic brane- world (without 

brane tension; U ~ models, first discussed by Binetruy, Deffayet and Langlois |Q, in which a(t) ~ t~, (which 
is valid at very high energies (p ^ A) as the initial singularity is approached; t ^ 0). Note that these solutions 
are self-similar, and are referred to as Friedmann brane- worlds Binetruy, Deffayet and Langlois solutions |^ or 
Brane- Robertson- Walker models |l^] . 

For the equilibrium point J^,, the 5 eigenvalues are: 

3(7 - 1), 3(7 - 1), (37 - 1), (37 - 1), (37 - 1) (40) 

Hence, for all physically relevant values of 7 (7 > 1), Th is a source (or past-attractor) in the brane-world scenario 
and the singularity is isotropic. This contrasts to the situation in GR in which anisotropy dominates for 7 < 2. This 
is also consistent with previous analyses of Bianchi type I and V models where Th is always a source for 7 > 1 (in the 
FRW models, Th is a source when 7 > 1/3 when U = ^ and 7 > 2/3 when ^ 0) ||l7[| . 



C. Discussion 



There are no other equilibrium points of the dynamical system that correspond to sources. However, it is instructive 
to consider the Kasner equilibrium points. 

The one-parameter set (circle) of (Bianchi I, shearing) Kasner vacuum (fi^ — 0) equilibrium points /C are given by 
Tj\ + Y?_ = \, and TVi = 7V2 = ^3 = 0. The eigenvalues are 

2(1 - 2E+), 2(1 + E+ + V3I]_), 2(1 + S+ - V3S_), -6(7 - 1), (41) 

and a fifth eigenvalue which is zero due to the fact that /C is a one-parameter set of equilibrium points. All of these 
equilibrium points are saddles; in particular, they can never be sources. The fourth eigenvalue is negative, whereas 
in GR it is positive, so that the structure of the stable and unstable manifolds close to K. is altered. In GR it was 
argued 0| that and A, defined by 

^={N^N2f ^{N^N^f ^{N^N{f, (42) 

satisfy il— >0, A^OasT— » —00, and hence in general all orbits approach K, (which is characterized by = 0, 
A = 0) into the past. A qualitative and numerical investigation indicates that does not tend to zero as r — > — 00, 
and hence generically K is not, nor does it constitute part of, a (past-) attracting set in the models under investigation. 

Close to the singularity we expect that 7 > 4/3. However, we note that there is a bifurcation at 7 = 1, and T\, and 
/C coalesce to form a two-parameter set of equilibrium points J which are analogues of the Jacobs stiff fluid solution 
in GR in which fif, 7^ 0. A subset of these equilibrium points are sources. In general -I- Y?_ ^ (the equilibrium 
points with zero shear are the special points in J equivalent to J-b)\ however, there are sources in J in which + 1]^ 
(and hence the shear) is arbitrarily small. 
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IV. CONCLUSIONS 



Therefore, generically the Bianchi IX branc-world models do not have space-Uke and oscillatory singularities in the 
past, and consequently brane- world cosmological models do not exhibit Mixmaster and chaotic- like behaviour close 
to the initial singularity. 

We expect this to be a generic feature of more general cosmological models in the brane-world scenario. In particular, 
we anticipate that !Fb is a source in all spatially homogeneous models. Following arguments similar to those of section 
III. A, it follows that there exists a singularity of a similar nature in all orthogonal Bianchi brane-world models. In 
general, fif, will again dominate as the initial singularity is approached into the past and the qualitative results of 
section III.B will follow. In particular, Th will be a local source and in general the initial singularity will be isotropic. 

A. Inhomogeneous Brane- World Models 

Let us consider the dynamics of a class of spatially inhomogeneous cosmological models with one spatial degree of 
freedom in the brane-world scenario. The G2 cosmological models admit a 2-parameter Abelian isometry group acting 
transitively on spacelike 2-surfaces. These models admit one degree of freedom as regards spatial inhomogeneity, 
and the resulting governing system of evolution equations constitute a system of autonomous partial differential 
equations in two independent variables. We follow the formalism of pO[ | which utilizes area expansion normalized 
scale-invariant dependent variables, and we use the timelike area gauge to discuss the asymptotic evolution of the 
class of orthogonally transitive G2 cosmologies near the cosmological initial singularity. The initial singularity can be 
shown to be characterized by .Ei^ — > as r — > — 00, where Ei^ is a normalized frame variable [2C[| . 

The 'equilibrium point' Tb is characterized by zero shear = Sx =0) and zero curvature (A^_ = A^x = 0) 
(where we note that the area expansion normalized variables utilized in [po| are equivalent to H-normalized variables 
close to Tb)- Linearizing the evolution equations about Tb^ using the same spatial reparametrisation as in [ pO[ (so 
that Ei^ — exp((37 — 1 )■'')), we obtain the following general solution of the linearized equations in which the shear 
variables satisfy 

E_ =ai(x)exp(3(7-l)T), Ex = a2(x) exp (3(7 - 1)t), (43) 
the curvature variables satisfy 

N. = az{x) exp ((37 - 1)t), A^x = a^{x) exp ((37 - 1)t), (44) 

and 

rib = 1 + 05(2;) exp(3(7- 1)t), (45) 

where the ai{x) are arbitrary functions of the space coordinate. As in |^^, it can be shown that the tilt or peculiar 
velocity v between the timelike area gauge and the fluid 4- velocity obeys v — a(,{x) exp ((37 — 1)t) and hence tends 
to zero as T — !■ —00 and, in addition, it follows from the conservation laws that J7 — s- 0, ilu ^ and JIa ^ as the 
initial singularity is approached. 

Strictly speaking, an inhomogeneous energy density is not consistent with — 0, V^v — 0; in particular, a 
non-zero D^p acts as a source for Q^. (Physically V^u corresponds to graviational waves and will not affect the the 
dynamics close to the singularity.) Writing | |= Q, an analysis of the evolution equation for Q close to Tb 
then yields QH'"^ ~ Q{x) exp (2(37 — 2)t), so that QH ^ — s- as r ^ —00 (as the initial singularity is approached). 
Moreover, it can then be shown that a small Q does not affect the dynamical evolution of U (to lowest order) close to 
!Fb (indeed, Q/U ~ 0(1), so that if Q is negligible it remains negligible). Hence, equations (43) - (45) are consistent 
with the evolution of Q^, and consequently represent a self-consistent solution close to Tb- 

The above linearized solution represents a general solution in the neighbourhood of the initial singularity. Hence 
Tb is a local source or past-attractor in this family of spatially inhomogeneous cosmological models for 7 > 1. In 
particular, we see that the shear and curvature asymptote to zero as r — > —00, and hence the singularity is isotropic. 
We also note that, unlike the analysis of the perfect fluid GR models in the Kasner equilibrium set K. are found 
to be saddles in the class of G2 brane-world cosmological models. 

The most detailed proposal for the structure of space-time singularities in GR are the conjectures of Belinskii, 
Khalatnikov and Lifshitz (BKL) one aspect of which is that each spatial point evolves towards the singularity as 
if it were a spatially homogeneous (Bianchi) cosmology. That, is, generic space-times have the property that spatial 
points decouple near the singularity, the EFE effectively reduce to ordinary differential equations (i.e., the spatial 



7 



derivatives have a negligible effect on the dynamics) and the local dynamical behaviour is asymptotically like that 
of Bianchi models near the singularity. Support for this conjecture comes from recent analyses that show that the 
presence of the inhomogeneity ceases to govern the dynamics asymptotically toward the singularity in particular classes 
of inhomogeneous models. In a special class of Abelian G2 spatially inhomogeneous models, the so-called "velocity- 
dominated" spacetimes, the evolution at different spatial points will approach that of different Kasner solutions [ p2[ . 
A numerical investigation of a class of vacuum Gowdy G2 cosmological spacetimes that represent an inhomogeneous 
generalization of Bianchi type VIq models with a magnetic field has shown evidence that at a generic point in space the 
evolution toward the initial singularity is asymptotically that of a spatially homogeneous spacetime with Mixmaster 
behaviour |Q . More rigorous results are given in |^,^ . 

The dynamics of the class of inhomogeneous brane-world models considered above, together with the BKL conjec- 
tures, indicate that a wide class of inhomogeneous brane-world models will have an isotropic initial singularity. 



B. Discussion 



Therefore, it is plausible that generically the brane-world cosmological models have an isotropic initial singularity, 
whose generic evolution near the cosmological initial singularity is approximated by a spatially homogeneous and 
isotropic model in a rigorously defined mathematical sense |p5f| (in contrast to the situation in GR). Therefore, it 
is possible that brane cosmology provides for a quiescent Eq] initial cosmological period in which the Universe is 
smooth and highly symmetric, perhaps explained by entropy arguments and the second law of thermodynamics [ p7| , 
thereby naturally providing the precise conditions for inflation to subsequently take place and consequently avoiding 
the problem of initial conditions in inflation [Ts| ] . 

The results in this paper are incomplete in that a description of the gravitational field in the bulk is not provided. 
Unfortunately, the evolution of the anisotropic stress part is not determined on the brane. These nonlocal terms also 
enter into crucial dynamical equations, such as the Raychaudhuri equation and the shear propagation equation, and 
can lead to important changes from the GR case. The correction terms must be consistently derived from the higher- 
dimensional equations. Additional modifications can also occur for higher-dimensions (than 5), for more general 
(than static) higher-dimensional (bulk) geometry, for higher (non-GR) curvature corrections, for higher-dimensional 
matter fields (e.g., scalar fields) in the bulk and for motion of the brane j28|. Further work is therefore necessary. 
However, it is plausible that the main results of this paper will persist when these additional effects are included. In 
particular, this is supported by a dynamical analysis of the evolution of and the fact that since Vi_ii, corresponds 
to gravitational waves in higher-dimensions it is expected that the dynamics will not be affected significantly at early 
times close to the singularity. 

It is also of importance to examine what happens when anisotropic stresses are included. Recently, Barrow and 
Hervik studied a class of Bianchi type I brane-world models with a pure magnetic field and a perfect fluid with 
a linear barotropic 7-law equation of state. They found that when 7 > |, the equilibrium point is again a local 
source (past-attractor), but that there exists a second equilibrium point denoted PHi, which corresponds to a new 
brane-world solution with a non-trivial magnetic field, which is also a local source. When 7 < |, PHi is the only 
local source (the equilibrium point !Fb is unstable to magnetic field perturbations and hence is a saddle); however, 
even in this case it was argued that chaotic behaviour is not possible ]30[ | . Finally, we note that Bianchi type I and IX 
models have also been studied in Hofava-Witten cosmology in which the fifth coordinate is a /Z2 orbifold while the 
remaining six dimensions have already been compactified on a Calabi-Yau space, and it was argued [ p9[ that there is 
no chaotic behaviour in such Bianchi IX Hofava-Witten cosmologies. 
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